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Director.  Dr.  J.  C.  Liebinan  is  Head  of  the  Department  of  Civil  Engineering  at  UIUC.  and 
Dr.  R.  Quattrone  is  Chief,  CERL-EM. 


CONTENTS 


DD  FORM  1473 
FOREWORD 

LIST  OF  TABLES  AND  FIGURES 

1  INTRODUCTION . 

Background 

Objective 

Approach 

2  METHODOLOGY  . 

Cap-Member-Soil  Model 
Cap-Member  Model 

3  CAP  MEMBER-SOIL  MODELS . 

Introduction 
Method  of  Analysis 
Soil  and  Concrete  Parameters 
Parametric  Analysis 

4  CAP-MEMBER  MODELS . 

Introduction 

Exploratory  Finite  Element  Models 
HP  14X117  Typical  Design  Loading  Analyses 
HP  14X117  Parametric  Analysis 

5  CONCLUSIONS . 

REFERENCES 


DISTRIBUTION 


TABLES 


Number 

1  Model  Parameters  and  Ranges 

2  Typical  Design  Loadings  for  HP  14X73  and  HP  14X117  Piles 

3  Values  of  nh  for  Sands 

4  So>l  Coefficients  for  Clays 

5  Limit  Elastic  Deformation  (y,e)  for  Sands 

6  Limit  Elastic  Deformation  (yle)  for  Clays 

7  Comparison  Between  Ihs  and  Irb 

8  Computed  Response  to  Typical  Design  Situations 

9  Angular  Rotation  of  HP  1 2X74  Pile 

10  Index  Correlating  Analysis  Case  and  Model  Parameters 

1 1  Basic  Parameters  for  HP  14X1 17  Parametric  Analysis 

FIGURES 

1  Model  Parameters  and  Ranges 

2  Design  Loadings  for  HP  14X73  and  HP  14X1 17  Piles 

3  Graphical  Definition  of  p  and  y 

4  Typical  Load-Deformation  Curves  (p-y  curves)  for  Sands  and  Clays 

5  Elastic  Half-Space  Under  a  Strip  Load.  Generalized  Solution 

6  Elastic  Half-Space  Under  a  Strip-Load,  Centralized  Solution 

7  Plot  of  IRB  in  Terms  of  m  =  4a  and  v  =  0.2 

8  Typical  Moment  and  Deflection  Curves  for  HP  14X1 17  Pile 
with  12  in.  Embedment 

9  Typical  Moment  and  Deflection  Curves  for  HP  14X117  Pile 
with  48  in.  Embedment 

10  Effect  of  Le./B  Ratio  on  Pile  Rotation  at  Soil-Cap  Interface, 

Typical  Design  Situation  for  HP  14X73  Pile 


5 


FIGURES  (Cont'd) 


Number 

1 1  Effect  of  Le/B  Ratio  on  Pile  Rotation  at  Soil-Cap  Interface, 
Typical  Design  Situation  for  HP  14X117  Pile 

12  Effect  of  Le/B  Ratio  on  Pile  Rotation  at  Soil-Cap  Interface, 
HP  14X73  Pile  with  15  kip  Shear  at  Base  of  Cap 

13  Effect  of  Le/B  Ratio  on  Pile  Rotation  at  Soil-Cap  Interface, 
HP  14X1 17  Pile  with  15  kip  Shear  at  Base  of  Cap 

14  Effect  of  Le/B  Ratio  on  Pile  Rotation  at  Soil-Cap  Interface, 
HP  14X1 17  Pile  in  Loose  Sand,  nh  =  4  ton/cu  ft 

15  Effect  of  Le/B  Ratio  on  Pile  Rotation  at  Soil-Cap  Interface, 
HP  14X117  Pile  in  Dense  Sand,  nh  =  34  ton/cu  ft 

16  Effect  of  Le/B  Ratio  on  Pile  Rotation  at  Soil-Cap  Interface, 
HP  14X73  Pile  in  Stiff  Clay,  k  =  50  tons/cu  ft 

1 7  Effect  of  Le/B  Ratio  on  Pile  Rotation  at  Soil-Cap  Interface, 
HP  14X73  Pile  in  Hard  Clay,  k  =  200  tons/cu  ft 

18  Effect  of  Le/B  Ratio  on  Pile  Rotation  at  Soil-Cap  Interface, 
Various  Pile  Sections  in  Medium  Dense  Sand  with 

5  kip  Shear  at  Base  of  Cap 

19  Effect  of  Le/B  Ratio  on  Pile  Rotation  at  Soil-Cap  Interface, 
Various  Pile  Sections  in  Medium  Dense  Sand  with 

10  kip  Shear  at  Base  of  Cap 

20  Stress  Distribution  in  an  HP  Pile  Section 

21  Exploratory  Finite  Element  Models 

22  Lateral  Displacement  Curves  for  HP  1 2X74  Pile  with 
0.5D,  D,  2D,  and  4D  Embedments 

23  Finite  Element  Representation  of  HP  Pile 

24  Two-Dimensional  Interface  Element 

25  Force-deflection  Relationship  for  the  Interface  Element 

26  Typical  Displacement  Field  for  HP  14X117  Pile 

27  Typical  Stress  Contour  Plot  for  ax 


28  Typical  Stress  Contour  Plot  for  o, 


FIGURES  (Cont'd) 


Number 

29  Typical  Stress  Contour  Plot  for  rxy 

30  Typical  Stress  Contour  Plot  for  amax 

31  Typical  Stress  Contour  Plot  for  omin 

32  Typical  Stress  Contour  Plot  for  rmax 

33  Lateral  Displacement  Curves  for  Centerline  of  HP  14X1 17  Pile 
with  2D  and  4D  Embedments,  /a  =  0.4 

34  Interface  Element  Forces  for  HP  14X117  Pile  with 
2D  Embedment,  /a  =  0.4 

35  Interface  Element  Forces  for  HP  14X1 17  Pile  with 
4D  Embedment,  n  =  0.4 

36  Lateral  Displacement  Curves  for  Centerline  of  HP  1 4X 1 1 7  Pile 
with  2D  Embedment,  n  -  0.4,  0.7,  and  1 .4 

37  Interface  Element  Forces  for  HP  14X117  Pile  with 
2D  Embedment,  /a  =  0.7 

38  Interface  Element  Forces  for  HP  14X1 17  Pile  with 
2D  Embedment,  /a  =  1 .4 

39  Sections  in  Concrete  Cap  Used  for  Normal  and  Shear  Force 
Distribution  Comparisons 

40  Normal  Force  Distribution  at  Section  A-A  for  Cap  Thickness  of  4D 

41  Shear  Force  Distribution  at  Section  A-A  for  Cap  Thickness  of  4D 

42  Normal  Force  Distribution  at  Section  B-B  for  Cap  Thickness  of  4D 

43  Shear  Force  Distribution  at  Section  B-B  for  Cap  Thickness  of  4D 

44  Lateral  Displacement  Curves  for  Centerline  of  HP  14X117  Pile  for 
Effective  Widths  of  Concrete  Cap  of  14.9  in.  and  44.7  in. 

45  Lateral  Displacement  Curves  for  Centerline  of  HP  14X117  Pile 
for  Pile  Spacings  of  3D  to  5D 

46  Normal  Force  Distribution  at  Section  A-A  for  Cap  Thicknesses 
of  4D  to  7D 


FIGURES  (Cont'd) 


Number  Page 

47  Shear  Force  Distribution  at  Section  A-A  for  Cap  Thicknesses 

of  4D  to  7D  59 

48  Normal  Force  Distribution  at  Section  B-B  for  Cap  Thicknesses 

of  4D  to  7D  60 

49  Shear  Force  Distribution  at  Section  B-B  for  Cap  Thicknesses 

of  4D  to  7D  61 


FIXITY  OF  MEMBERS  EMBEDDED  IN 
CONCRETE 


1  INTRODUCTION 


Background 

The  Corps  of  Engineers  designs  and  constructs  many 
concrete  structures  with  pile  foundations  such  as  locks, 
dams,  pumping  stations,  bridges,  and  lloodwalls. 
Structural  piles  (steel,  concrete,  or  timber)  are  used  to 
transfer  the  horizontal  and  vertical  loadings  on  these 
structures  to  soil  layers  which  can  provide  adequate 
support.  The  degree  of  fixity  of  the  piles  is  an  im¬ 
portant  assumption  in  the  analysis  and  design  of  the 
pile  groups  that  support  the  structures.  For  instance, 
elastic  solutions  for  a  pile  subjected  to  lateral  load  and 
moment  forces  indicate  that  a  structure  loaded  hori¬ 
zontally  will  deflect  approximately  one-quarter  to  one- 
half  as  much  for  piles  that  are  fully  fixed  at  the  base  of 
the  structure  as  opposed  to  piles  that  are  simply  pinned 
to  the  base  of  the  structure.1  For  concrete  structures 
founded  on  piles,  the  resistance  to  lateral  forces 
depends  on  the  degree  of  fixity  of  the  embedded  piles 
and  the  complex  interaction  of  the  piles  with  the 
surrounding  soil.  For  example,  piles  embedded  a 
limited  depth  will  resist  only  shear  and  axial  loads, 
piles  embedded  an  adequate  depth  will  resist  moments, 
shears,  and  axial  loads  and  significantly  reduce  the 
lateral  deflections  for  the  siructure. 

Pile  foundations  represent  a  considerable  portion  of 
the  total  cost  of  a  structure  founded  on  piles.  Conse¬ 
quently.  the  length  of  pile  embedment  required  to  de¬ 
velop  a  condition  of  full  fixity  has  a  significant  impact 
on  the  cost  of  the  pile  foundation  and  ultimately  the 
cost  of  the  structure.  Corps-wide,  this  cost  could  be  in 
the  millions  of  dollars. 

A  “rule  of  thumb"  is  commonly  used  in  the  design 
for  steel  HP  shape  piles  for  determining  embedment 
length.  An  embedment  of  1  ft  or  less  into  the  concrete 
structure  is  considered  pinned,  a  1-  to  2-ft  embedment 
is  considered  partially  fixed,  and  a  4-ft  embedment  is 
considered  fully  fixed.  However,  there  are  no  known 
theoretical  or  experimental  data  that  support  these  as¬ 
sumptions.  For  some  structures,  it  is  desirable  to  fully 


lH.  (I.  Poulos  and  I  .  H.  Davis.  Mastic  Solutions  for  Soil 
and  Rock  Mechanics  (John  Wiley  and  Sons.  Inc..  1974). 


fix  the  piles  to  limit  deflection.  A  considerable  amount 
of  theoretical  and  experimental  research  is  available  to 
define  the  load  transfer  mechanism  tor  a  pile  em¬ 
bedded  in  the  soil  below  the  pile  cap;  however,  there 
appears  to  be  little  data  available  to  define  the  load 
transfer  mechanism  for  the  portion  embedded  in  the 
concrete  pile  cap. 

At  the  request  of  the  St.  Louts  District.  US  Army 
Corps  of  Engineers,  a  two-part  study  was  performed  to 
develop  relationships  between  the  degree  of  fixity  and 
the  pile  embedment  length.  The  studies  were  per¬ 
formed  through  a  joint  effort  of  the  U.S.  Army  Con¬ 
struction  Engineering  Re.  :  Laboratory  (CERL) 
and  the  University  of  II'  s  at  Urbana-Champaign 
(UIUC).  Department  of  Cb  ''’ineering. 

Objective 

The  objective  of  the  ,  ,r  research  is  to  deter¬ 
mine  the  pile  embedment  ..  .gth  required  to  provide 
reasonable  assurance  that  a  condition  of  full  fixity 
exists  for  the  pile  embedded  in  the  base  of  the 
structure  and  that  acceptable  stresses  created  by  the 
moment  developed  in  the  restrained  condition  are 
transferred  from  the  pile  to  the  concrete  cap  and  vice 
versa 

Approach 

CTRL  and  UIUC  examined  the  fixity  of  embedded 
piles  through  the  use  of  computer  models.  UIUC  used 
cap-mcmbcr-soil  models  in  its  analysis;  CERL  ex¬ 
amined  the  structural  integrity  of  piles  using  cap- 
member  models.  Parameters  and  ranges  used  in  these 
studies  were  developed  by  CERL  and  coordinated  with 
the  St.  Louis  District.  A  more  detailed  explanation  of 
the  approach  used  in  these  studies  is  provided  in  the 
following  chapter  on  methodology. 


2  METHODOLOGY 


The  problem  of  what  pile  length  will  ensure  full 
fixity  is  analogous  to  a  beam-column  connection  which 
is  usually  assumed  to  be  perfectly  rigid.  But  if  a  more 
precise  determinativ-n  of  the  deflections  is  made,  con¬ 
sidering  deformation  of  the  joint,  the  angles  existing 
between  connecting  members  at  the  joints  will  not  be 
preserved.  Likewise  for  steel  piles  embedded  in  a  con¬ 
crete  cap:  the  angle  between  the  pile  and  concrete  cap 
will  not  be  maintained  after  loading  has  been  applied 
to  the  structure  due  to  the  deformability  and  degrada¬ 
tion  of  the  materials  involved,  especially  the  concrete. 


Analytical  studies  were  performed  by  CFRL  and 
UIUC  using  models  of  the  cap-member  system  and  the 
cap-member-soil  system.  These  studies  were  initially 
conducted  using  tire  parameters  and  ranges  presented 
in  Table  I  and  diagrammed  in  Figure  1.  These  parame¬ 
ters  and  ranges  were  developed  by  CERL  in  conjunc¬ 
tion  witit  the  St.  Louis  District  before  this  investigation 
began.  Preliminary  results  from  both  the  CERL  and  the 
UIUC  studies  using  these  parameters  indicated,  how¬ 
ever.  that  an  embedment  length  of  approximately 
twice  the  diameter  or  depth  of  the  pile  was  adequate  to 
develop  a  condition  approximating  full  fixity.  The  St. 
Louis  District,  therefore,  requested  CERL  and  UIUC  to 
investigate  the  behavior  of  IIP  14X1  17  and  HP  14X73 
steel  piles  with  2-ft  and  4-ft  embedments  in  a  concrete 
cap  and  subjected  to  the  typical  design  loadings  condi¬ 
tions  presented  in  Table  2  and  Figure  2. 


Table  1 

Model  Parameters  and  Ranges 

Parameter  Range 

Member  Spacing  3  to  5D* 

Member  I  mbedment  1  to  4D 

Cap  Thickness  6  to  10D 

Member  Configuration  I D I 
Circular  In  to  24  in  in  dia 

“H”  in  .  10  to  14  /  14': 

Modulus  ol  I  last ici t \ 

Concrete 
Steel 

Poisson’s  Ratio 

Concrete  0.15  to  0.20 

Steel  0.3 

Soil  C  1  to  0.5 


3  to  5  /  1 0'  psi 
29.6  >  10'  psi 


Cap-Member-Soil  Model 

UIUC  represented  the  concrete  cap-pile-soil  system 
as  a  Winkler  model,  i.e.,  the  pile  was  supported  by  a 


Soil  Subgrade  Reaction 
Sands 
Clays 


4  to  56  tons/cu  It 
50  to  200  tons/cu  ft 


series  of  independent  springs.2  Springs  with  linear  load-  Soil  I  mbedment  Length  20  to  60  ft 

deformation  characteristics  represented  the  concrete 

cap,  and  springs  with  nonlinear  characteristics  repre-  — - - - 

sented  the  soil.  The  pile  was  assumed  to  be  linearly  *D  designates  member  diameter  or  depth  depending  on  niettt- 
........  ..  ,  .  her  configuration. 

elastic.  Variations  ot  the  significant  parameters  were  in¬ 
vestigated  to  study  the  relationships  among  tire  degree 

of  fixity  of  the  pile,  the  pile  embedment  length  in  the 
3L  C.  Reese.  ’‘Laterally  Loaded  Piles:  Program  Documen-  concrete  cap.  and  the  interaction  of  the  pile  with  the 

tation."  ASCI:  PnH'vedinns.  Vol.  103.  No.  CT4  (April  1977),  surrounding  soil.  UIUC  used  the  computer  program 

pp.  287-305  COM  622  to  perform  the  calculations  for  the  cap- 
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Figure  1 .  Model  parameters  and  ranges. 


Table  2 

Typical  Design  Loadings  for  HP  14X73  and  HP  14X1 17  Piles 


Y, 

K4 

F 

E, 

(kips) 

(kips) 

(kips) 

(in. -kips) 

(in. -kips) 

(in. -kips) 

NORMAL  OPKRATINC. 

Pinned  HP  14X73  nj,  ~  .017  kips/eu  in 

13 

6.4 

125 

0 

0 

0 

Lnd  HP  14X117  nh  -  .003 

0 

0 

188 

0 

0 

II 

fixed  HP  14X73  nh  =  .001 

2.8 

0 

191 

0 

295 

II 

1-nd  HP  14X117  nh  =  .001 

2.2 

0 

138 

0 

260 

0 

NORMAL  OPERATING 
<0.1  g  E.Q.) 

Pinned  HP  14X73  nh  =  .027 

-26 

-1.7 

123 

0 

0 

0 

End  HP  14X117  nh  =  .03 

22 

1 

201 

n 

0 

0 

fixed  HP  14X73  nh  =  .03 

18 

0 

228 

0 

954 

0 

End  HP  14X117  nh  =  .03 

15 

0 

142 

0 

913 

0 

member-soil  model.3  The  program  was  modified  to  To  minimize  computer  costs  and  to  produce  results 

model  the  embedment  of  the  pile  in  the  concrete  cap.  that  were  readily  comparable  with  the  L'lL'C  results. 

To  accurately  represent  the  different  soil  conditions  it  only  two-dimensional  linear  and  nonlinear  models  of 

was  necessary  to  determine  soil  spring  constants  for 
sands  and  clays  of  various  strengths  and  limit  pressures. 

Spring  constants  were  also  determined  for  the  concrete 
cap  section  of  the  model. 

The  relationships  among  the  parameters  were  analy¬ 
zed  for  various  values  of  soil  stiffness,  pile  embedment 
length,  pile  type,  and  shear  force  at  the  interface  be¬ 
tween  the  soil  and  the  cap. 

Cap-Member  Model 

CERL  analyzed  the  problem  from  a  purely  struc¬ 
tural  standpoint.  Finite  element  models  of  the  pile 
embedded  in  the  concrete  cap  were  assembled  without 
consideration  of  the  underlying  soil  material.  Varia¬ 
tions  of  the  significant  parameters  were  investigated  in 
an  effort  to  determine  guidance  on  the  length  of  pile 
embedment  required  to  produce  a  condition  approxi¬ 
mating  full  fixity. 

A  general  purpose  computer  program,  ANSYS,  was 
used  to  analyze  the  embedment-fixity  problem. 

ANSYS  is  a  large-scale  finite  elements  program  that  can 
analyze  two-dimensional  and  three-dimensional  linear 
or  nonlinear  models  under  static  and  dynamic  condi¬ 
tions,  and  small  or  large  displacements.4 

3Reese.  pp  287-305. 

A  ANSYS  User's  Manual  Revision  3  (Swanson  Analysis  Sys-  Figure  2.  Design  loadings  for  HP  14X73  and  HP 
terns.  Inc.,  Houston.  Pennsylvania,  1 977).  14X1  17  piles. 
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the  cap-member  system  were  analyzed.  Dimensions  for 
the  two-dimensional  models  were  based  'in  the  ranges 
in  Table  I  and  equal  orthogonal  spacing  01  the  mem¬ 
bers  was  assumed.  The  finite  element  used  in  the  two- 
dimensional  linear  models  was  STIF  42.  a  general 
purpose  element  for  modeling  solid  structures.  The 
element  was  used  in  a  biaxial  plane  strain  state.  The 
element  is  defined  by  four  nodal  points,  each  having 
two  degrees  of  freedom  translations  in  the  nodal  x 
and  y  directions.  For  the  two-dimensional  nonlinear 
models,  the  STIF  42  element  was  also  used;  however, 
the  flanges  of  the  pile  were  modeled  with  a  plane  bar 
element.  STIF  1.  The  interface  between  the  pile  sec¬ 
tion  and  the  concrete  cap  was  modeled  with  a  two- 
dimensional  interface  element.  STIF  12.  The  interface 
element  represents  two  plane  surfaces  which  may  main¬ 
tain  or  break  physical  contact  and  which  may  slide 
relative  to  each  other.  The  element  is  capable  of  sup¬ 
porting  only  compressive  forces  in  the  direction  nor¬ 
mal  to  its  surfaces,  and  shear  (Coulomb  friction)  in  its 
tangential  direction.  The  element  has  two  degrees  of 
freedom  at  each  node  translations  in  the  nodal  x  and 
y  directions.  The  input  properties  required  to  repre¬ 
sent  the  model's  materials  characteristics  were  the 
modulus  of  elasticity.  Poisson’s  ratio,  coefficient  of 
friction,  and  shear  modulus. 

The  measure  of  pile  fixity  is  based  on  rotational 
response  of  the  embedded  pile  section  at  the  concrete 
cap  face.  Since  the  ANSYS  STIF  42  element  provides 
only  displacements  at  its  nodal  points,  a  procedure  was 
developed  to  interrogate  the  ANSYS  displacement  re¬ 
sults  and  obtain  displacement  values  for  selected  nodes. 
These  displacements  were  then  used  to  compute  the 
required  rotational  values  at  the  cap  face  for  the 
various  cases  analyzed. 


3  CAP  MEMBER-SOIL  MODELS 


Introduction 

In  the  UIUC  analytical  study,  the  degree  of  fixity  of 
the  pile  embedded  in  a  concrete  cap  and  the  underlying 
soil  was  analyzed  as  two  independent  pile  sections,  i.e., 
one  in  the  cap  and  the  other  in  the  soil.  The  concrete 
cap  was  modeled  as  springs  with  linear  load-deforma¬ 
tion  characteristics,  and  the  soil  was  represented  by 
springs  with  elastic-plastic  load-deformation  charac¬ 
teristics.  The  pile  was  assumed  to  be  linearly  clastic. 
UIUC  modified  the  COM  622  computer  program  to 
perform  the  analysis. 


Method  of  Analysis 

Original  Method 

The  COM  622  program  was  modified  to  solve  the 
problem  of  a  beam  resting  on  a  nonlinear  Winkler's 
foundation  with  a  finite  difference  numerical  scheme 
The  boundary  conditions  (four  conditions  are  required 
to  solve  the  problem)  are  preprogrammed  at  the  tip  ol 
the  pile  and  are  input-dependent  at  the  head  of  the  pile 
(soil-cap  interface). 

The  conditions  of  zero  moment  and  zero  shear  at 
the  tip  of  the  pile  are  preprogrammed  in  COM  622.  At 
the  head  of  the  pile,  the  input  boundary  conditions  can 
be  shear  and  moment,  shear  and  slope,  or  shear  and 
ratio  of  moment  to  slope  (rotational  spring  constant ). 

Modified  Method 

TIk  approach  taken  by  UIUC  consists  of  a  two-step 
iterative  scheme  which  was  used  because  only  one  of 
the  two  boundary  conditions  at  the  head  of  the  pile 
was  known.  The  full  pile  is  modeled  in  the  cap  and  in 
the  soil  as  two  independent  pile  sections.  The  analysis 
is  carried  out  for  each  section  separately.  The  head  of 
each  pile  section  begins  at  the  soil-cap  interface.  Itera¬ 
tions  arc  performed  until  the  continuity  of  boundary 
conditions  at  the  soil-cap  interface  is  achieved. 

A  complete  iterative  cycle  consists  of  the  following 
steps: 

Step  one;  Analysis  of  pile  section  in  concrete  cap. 

(a)  Assign  moment:  the  moment  is  assumed  to  be 
zero  at  the  head  of  the  pile  for  the  first  iteration.  For 
the  second  iteration,  the  moment  assigned  to  the  head 
is  the  value  computed  in  (d)  below  for  the  first  itera¬ 
tion.  Note  that  the  shear  in  the  pile  at  the  concrete  soil 
interface  is  a  known  value  for  the  problem. 

(b)  Perform  analysis:  compute  the  slope  of  the  cap- 
pile  section  at  the  soil-cap  interface  for  the  value  of  the 
assigned  moment  and  shear  from  (a). 

Step  two:  Analysis  of  pile  section  in  soil. 

(c)  Assign  the  slope  computed  in  (b). 

(d)  Perform  analysis:  obtain  the  computed  value  of 
moment  in  the  soil-pile  section  at  the  soil-cap  interface. 

(e)  Compare  the  moment  computed  in  (d)  with  the 
moment  assigned  in  (a).  If  the  two  values  aie  not 
within  an  acceptable  tolerance,  perform  a  new  iterative 


cycle.  Otherwise,  the  solution  has  been  reached. 
(Acceptable  tolerance  is  reached  when 


of  elasticity  of  the  sand  for  a  given  depth  /  can  be 
defined  as 


M,  -M2 

M, 


X  100 


<  1 


where  M,  and  \1:  are  moments  computed  from  cap- 
pile  and  soil-pile  analysis.) 


Es  =  X  av  =  X  7  / 

where  X  is  a  coefficient  which  depends  only  on  the 
density  of  the  sand,  and  ov  is  the  effective  overburden 
pressure  at  the  depth  /. 


The  iterative  procedure  is  automatic  in  the  modified 
version  of  the  COM  622  computer  program.  A  degree 
of  convergence  acceptable  for  engineering  purposes 
was  usually  achieved  within  four  cycles  of  iterations. 

Soil  and  Concrete  Parameters 

For  the  solution  to  the  problem  of  a  laterally  loaded 
pile  by  the  method  used  in  the  COM  622  computer 
program,  it  is  necessary  to  predict  a  set  of  load-defor¬ 
mation  curves  (known  as  p-y  curves)  for  each  medium 
(soil  and  concrete).  Figure  3  shows  the  graphic  defini¬ 
tion  of  p  and  y  and  a  typical  set  of  p-y  curves.  The 
following  two  sections  deal  with  construction  of  p-y 
curves  for  soil  and  concrete. 

Soil  Parameters 

Sands  and  clays  are  the  two  soil  types  considered  in 
this  study.  For  pressure  levels  below  the  limit  pressure, 
the  stiffness  values  used  are  from  Terzaghi.5  The 
method  described  below  was  used  to  establish  limit 
pressures  for  both  sand  and  day.  Klastic-plastic  p-y 
curves  were  constructed  for  both  sands  and  clays  by 
combining  Terzaghi's  recommendation  for  the  elastic 
portion  and  using  the  limit  pressure  calculations  to 
establish  the  upper  limit  of  p. 

The  derivations  and  values  retained  for  analysis  are 
shown  in  the  following  subsections.  The  soil  coeffi¬ 
cient  (spring  constant)  characterizes  the  elastic  part 
while  the  limit  pressure  describes  the  plastic  region. 

Evaluation  of  Soil  Coefficients.  Soil  coefficients 
were  evaluated  for  sands,  clays,  and  a  typical  design 
situation  using  Terzaghi's  recommendations. 

1 .  Sands  The  soil  coefficient  k  for  sands  at  a  given 
depth  /.  below  the  surface  depends  on  the  width  B  of 
the  pile  (measured  at  a  right  angle  to  the  direction  of 
the  displacement),  the  effective  unit  weight  y  of  the 
sand,  and  the  relative  density  of  the  sand.  The  modulus 


5K.  Terzaghi.  "Kvaluation  of  Coefficients  of  Subgrade 
Reaction.”  (leoteehniquc.  Vol.  5.  No.  4  ( 1955).  pp  297-326. 


The  soil  coefficient  for  a  vertical  beam  with  width  B 
is  determined  by  the  relation  between  the  contact  pres¬ 
sure  p  and  the  corresponding  lateral  displacement  y  of 
the  contact  face.  The  theory  of  elasticity  defines  the 
relation  between  p  and  y  as 

F-s  _  X  7  > 

P  =  V  1 .35  B  ”  V  1 .35  B 
Hence  the  soil  coefficient: 


and  the  load  per  unit  length  of  the  pile: 

P  =  p  B  =  nh  t  y 

where  the  factor 

X  7 


The  value  of  X  ranges  from  about  100  for  \er\  loose 
sand  and  2000  for  dense  sand.  The  average  unit  weight 
of  dry  or  moist  sand  is  considered  as  y  -  100  lb  cu  ft. 
and  the  average  effective  unit  weight  of  submerged 
sand  as  y  =  60  Ib/cu  ft.  The  values  of  nh  given  in 
Table  3  are  consistent  with  these  values  of  X.  y,  and  y'. 

2.  Clays -For  piles  embedded  in  clays,  the  value  of 
the  soil  coefficient  k  increases  approximately  in  simple 
proportion  to  the  unconfined  compressive  strength  ot 
the  ciay  (qu).  Table  4  shows  the  proposed  value  of  soil 
coefficients  for  piles  embedded  in  clay. 


Table  3 

Values  of  %  for  Sands  (tons/cu  ft) 


Sand 

Loose 

Medium 

Dense 

(dry  nr  moist) 

7 

21 

56 

nh  (submerged) 

4 

14 

34 

Ground  Surfoce 


Graphical  Definition  of  Pand  y 


(a)  Side  View 

(b)  A-A,  Earth  Pressure  Distribution 
Prior  to  Lateral  Loading 

(c)  A-A,  Earth  Pressure  Distribution 
After  Lateral  Loading 

(d)  Set  of  P-y  Curves 


Figure  3.  Graphical  definition  of  p  and  y. 


Table  4 

Soil  Coefficients  for  Clays 
Very 

Clay  Stiff  Stiff  Hard 

q„  (tnn/sq  III  1-2  2-4  >4 

k  (tim/cu  ID  50  100  200 


3.  Typical  design  situation  Under  normal  static 
operation  conditions,  the  soil  was  represented  by  a 
very  loose  sand  with  nh  =  0.001  kips/cu  in.  =  1  ton/cu 
ft.  Under  normal  dynamic  operating  conditions  (0  1  g 
earthquake  acceleration),  the  soil  was  represented  by  a 
medium  dense  sand  with  nh  =  0.03  kips/cu  in.  =  26 
tons/cu  ft. 

Evaluation  of  Limit  Pressures.  The  limit  lateral  pres¬ 
sure  on  a  pile  undergoing  lateral  displacement  is  equal 
to  the  difference  between  passive  and  active  pressure. 
Sandy  soils  were  assumed  to  be  fully  submerged  (7'  = 
buoyant  unit  weight).  Clayey  soils  were  assumed  to 
respond  to  an  undraincd  loading  condition  (0  =  0 
analysis). 

1 .  Sands  Tire  limit  load  per  unit  length  of  pile  is 
given  by 


Table  5 


Limit  Elastic  Deformation  (y  le)  for  Sands 

Sandy 

Very 

Medium 

Soils 

Loose 

Loose 

Medium 

Dense 

Dense 

Typical 

Design 

Situation 

X 

X 

(tons/cu  ft) 

0.864 

4 

14 

26 

34 

<t>  (deg) 

30 

30 

37 

40 

45 

KP 

3.0 

3.0 

4.0 

4.6 

5.83 

kA 

0.33 

0.33 

0.25 

0.22 

0.17 

HP  14X117 
b  14.9  in. 

jy  (Ib/tn.) 

1.38 

1.38 

1.95 

2.27 

2.93 

>  ic  tin.) 

1.38 

0.298 

0.120 

0.076 

0.074 

HP  14X73 
b  =  14.6  in. 

j  (Ib/in.) 

1.35 

1.35 

1.91 

2.22 

2.87 

yie  fin.) 

1.35 

0.292 

0.1  18 

0.074 

0.07.3 

*z  measured  in  inches. 


P  =  (Kp  —  Ka  )y  /B 

where 

Kp  =  coefficient  of  passive  earth  pressure  ---  tan 
(45°  +0/2) 

KA  =  coefficient  of  active  earth  pressure  =  tan  (45° 

-  0/2) 

7'  =  buoyant  unit  weight  =  7^,  7^  =  (>0  lb 

cu  ft 

/  =  depth 
B  =  pile  width. 

The  ratio  p/z  =  (Kg  -  KA)7’B  is  a  characteristic  of 
both  soil  and  pile  width.  The  limit  elastic  deformation 
(yle)  is  computed  from  the  expression 

P  =  nh  /.  y 

which  yields: 

PI  ,1 

y,c  -  -  —  -  (Kp-KA)7B  — 

/.  nh  nh 

The  values  of  y  ,e  for  various  types  of  sands  and  two 
different  piles  (HP  14X73  and  HP  14X1  I7)arc  shown 
in  Table  5. 
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2.  Clays  In  an  undrained  loading  condition  (0  = 
0),  the  difference  between  passive  and  active  pressure 
(A)  is  equal  to  twice  the  uncontlned  compressive 
strength,  qu : 

A  =  2  qu 

The  limit  load  per  unit  length  of  pile  is  given  by 
P=  AB 

where  B  is  the  pile  width. 

The  limit  elastic  deflection.  y,c.  is  computed  from 
the  expression 

P=Kyu. 

or 

P 

yte  -  K 
or 


Table  6  shows  the  values  of  ylt  for  various  con¬ 
sistencies  of  clays  and  for  two  different  piles. 

Figure  4  shows  typical  load-deformation  curves  (p-y 
curves)  for  sands  and  clays. 

Determination  oj  Concrete  Spring  Constant 

Two  approaches  were  taken  to  evaluate  the  concrete 
spring  constant : 


Table  6 

Limit  Elastic  Deformation  (ylc)  for  Clays 


Clayey 

Soils 

Stiff 

Very 

Stiff 

Hard 

Uu  (lon/sq  it) 

1 

T 

3  ( low . 

k  (ton/cu  It) 

50 

100 

but  used  in 
unjlysis) 
200 

HP  14X117 

p  (lb/in.) 

414 

838 

1.342 

yic  <in-> 

0.597 

0.596 

0.447 

HP  14X73 

p (Ib/in.) 

406 

81  1 

1.217 

yte  (in  > 

0.585 

0.584 

0  438 

1.  Half-Space  Solution  The  equation  yielding  to 
compression  of  the  surface  layer  of  an  elastic  half¬ 
space  under  the  centerline  of  a  strip  load  was  derived. 
It  is  characterized  by  an  influence  coefficient.  IHS. 

2.  Elastic  Layer  over  Rigid  Base  Solution  The  nu¬ 
merical  solution  to  this  problem,  obtained  by  Poulos.6 
is  characterized  by  an  influence  coefficient,  1RB. 

It  was  found  that  both  coefficients.  lHs  and  IRB. 
are  approximately  equal  for  a  layer  thickness  at  least 
equal  to  a  width  and  a  half  of  the  loaded  area.  In  other 
words,  either  solution  can  be  used  in  evaluating  the 
concrete  spring  constant  as  long  as  the  representative 


f'H.  (j.  Poulos.  “Stresses  and  Displacements  in  an  Plastic 
Layer  Underlain  by  a  Rough  Rigid  Base,"  (ieolechniquc.  Vol. 
17,  No.  4  (1967),  pp  378-410. 


yfc  yic 


SAND 


CLAY 


Figure  4.  Typical  load-deformation  curves  (p-y  curves)  for  sand-  and  clays. 


concrete  thickness  on  either  side  of  a  pile  is  at  least 
equal  to  one  and  a  halftimes  the  width  of  the  pile.  The 
value  Kg  =  3084  kips/sq  in.  was  selected  for  all  the 
analyses.  It  was  computed  for  a  ratio  of  concrete  thick¬ 
ness  to  pile  width  equal  to  1.75  and  was  found  to  be 
approximately  equal  to  the  average  spring  constant  for 
values  of  concrete  thickness  to  pile  width  ratio  ranging 
from  one  to  five. 

Detailed  computations  of  K§  for  both  approaches 
are  given  in  the  following  pages. 

Half-Space  Solution.  Stresses  in  an  elastic  half-space 
under  a  strip  load  are:7 

a,  =  ^  [a  +  sina  cos(a  +  25)] 
a*  =  ■-  [a  +  sina  cos(a  +  26 )] 

-q 

a..  =  —  voc 

'  71 

q 

T\i  ~  ~  [sina  sin  (a  +  26)] 

where  ax,  ay,  and  a,  are  normal  stresses  in  x,  y,  and  /. 
directions  and  rx/  is  the  shear  stress  acting  in  /.  direc¬ 
tion  on  a  plane  normal  to  the  x  axis.  Other  variables  in¬ 
volved  in  the  above  equations  are  defined  in  Figure  5. 
Strain  in  /.  direction  can  be  defined  as: 

1 

f/  =  ^(0/  -  t'ffx  -  VO y) 


Figure  5.  Elastic  half-space  under  a  strip  load,  gen¬ 
eralized  solution. 

1  +  rHr 

e,  =  — - - a(  I  —  2n)  +  sina 

L  n 

From  Figure  6.  z  can  be  written  as 

a 

l  =  a  cot  - 

By  differentiation 

a  I 

d/.  =  -  - -  da 

-  .  2a 

sin^  - 

and  infinitesimal  settlement  dw  can  be  defined  as 


On  the  strip  centerline  (x  =  0)  ox,  ay,  a,.  rxz  simplify 
as  follows: 


_  Figure  6.  Elastic  half-space  under  a  strip  load,  ecn- 

7Poulosand  Davis  (1974).  tralized  solution 


v*  \vV  .vs.  %* 


lij 


dw  =  e,  d z 

Substituting  e,  and  dz  from  the  above  equations,  dw 
can  be  written  as 


q  I  +  c  — [a(l  —  2d)  +  sinal 

dw  =  - - —  a - da 

2n  L.  ,a 

sin  — 


Integration  of  dw  between  a  =  n  (z  =  0)  and  af 
2  Arc  cot  (m)  yields: 

rat  ,  q  1  +  v  , 

w=X  dw=^— 

X  [2a,  cot  ~  —  4  In  sin  ^-]  —  4  In  sin  ^-] 
replacing  a  by  L/ni  and  simplifying. 


2q  1  +  v  , 

w  = - - —  L[(  1  —  2d)  Arc  cot  (m) 

n  t 


+  (1  —  t>)  —  In  (1  +  m2)] 
m  v 


The  equivalent  spring  constant,  Ks,  is 


IRg  OS 


0.1  02  0.3  0L4  0.5 

2.0  15  1.0  05  0 

-L  m 

m 

Figure  7.  Plot  of  IRB  in  terms  of  m  =  4a  and  v  =  0.2. 


Table  7 

Comparison  Between 
IHs  and  Irb 


1  +  v  >/(l  2i')m  Arc  Cot  (m)  +  ( I  (v  ln(l  ■*  m2 )) 

The  settlement  of  tlte  surface  layer  of  an  elastic  half 
space  can  be  expressed  as 

2qL 

w  =  ^F*hs 


Ihs  =  ( 1  +  ^)[(  1  —  2d)  Arc  cot  (m ) 


+  (1  -n)-ln(l  +  m2)] 
m 


Elastic  Layer  over  Rigid  Base  Solution.  A  numerical 
solution  for  the  problem  of  the  settlement  of  a  finite 
layer  over  a  rigid  base  has  been  obtained  by  Poulos.8 
The  solution  can  be  expressed  as 


2qL 

w  =  ~^irb 


where  IRB  is  shown  in  Figure  7  in  terms  of  m  =  L/a  for 
a  Poisson's  ratio  of  0.2.  Table  7  shows  a  comparison 
between  IHS  and  1RB  for  various  m  =  L/a  values.  As 
can  be  seen  from  the  table,  both  methods  predict  very 
similar  settlements  of  the  layer  under  a  strip  load  for 
m  ^  2.0. 

Calculation  of  Concrete  Spring  Constant.  Using  the 
half-space  solution,  the  spring  constant  Ks  is  given  as 

_  2qa  _  2qL  rrE  1 

s  w  mw  m  Iuc 


’Poulos.  pp  378-410. 


18 


where 


and  concrete's  Young's  modulus,  from  ACI:9 

E(lb/st|  in)=  57000  VC  (lb/sc[  in) 

using  =  3000  lb/sq  in.,  then  E  =  3122  kips/sq  in.; 
and  concrete’s  Poisson’s  ratio 

v  =  0.2 

Thus,  concrete  spring  constant  Ks  with  m  =  3.5  is 
computed  as 

1hs=0.91 

9  ACI  318-77,  Building  Code  Requirements  for  Reinforeed 
Concrete  (American  Concrete  Institute  [ACI | ,  June  1978). 


2.  Soil  stiffness  and  soil  type 

3.  Shear  force  at  the  interface  between  the  soil  and 
cap 

4.  Pile  stiffness. 

Figures  8  and  9  show  a  typical  analysis  for  specific 
parameters.  The  analysis  was  performed  for  variation 
of  all  the  above  parameters;  the  results  are  shown  in 
Figures  10  through  19  and  Table  8. 

Figures  10  through  19  depict  the  rotation  of  pile  at 
the  soil-cap  interface  (0)  versus  the  ratio  of  embed¬ 
ment  length  in  the  cap  to  pile  width  (Le/B).  As  can  be 
seen  from  these  figures,  when  Le/B  is  less  than  two.  the 
rotation  of  the  pile  increases.  When  Le/B  is  more  than 
two,  the  rotation  of  the  pile  reaches  a  minimum 


Table  8 

Computed  Response  to  Typical  Design  Situations 

Fixed  End 


Embedded 

Shear  at 

Cap 

Computed 

Usual  Design 

Length 

Base  of  Cap 

Displacement 

Max.  Moment 

Moment 

Operation 

Soil 

Pile 

(in.) 

(kips) 

(in.) 

(in.-kips) 

(in.-kips) 

Normal 

Very  Loose 

HP  14X73 

12 

2.8 

0.27 

240 

295 

Static 

Sand 

24 

2.8 

0.23 

288 

36 

2.8 

0.23 

291 

48 

2.8 

0.23 

291 

nh  =  0.001  kip/cu  in. 

HP  14X1 17 

12 

2.2 

0.20 

202 

260 

^  1  ton/cu  ft 

24 

2.2 

0  16 

259 

36 

2.2 

0.15 

264 

48 

2.2 

0.15 

264 

Normal 

Medium 

HP  14X73 

12 

18 

0,44 

868 

954 

Dynamic 

Dense  Sand 

24 

18 

0.36 

1.140 

0.1  (! 

36 

18 

0.35 

1.160 

Earthquake 

48 

18 

0.35 

1.160 

nj,  =  0.03  kip/cu  in 

HP  14X117 

12 

15 

0.26 

605 

913 

~  26  ton/eu  ft 

24 

15 
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TYPICAL  DESIGN  SITUATION 
PILE  HP  14  x73 


EFFECT  OF  Le/B  RATIO  ON  PILE  ROTATION 
AT  SOIL-CAP  INTERFACE 


Rotio  of  Embedment  Length  to  Pile  Width,  Le/B  (in. /in.) 


Figure  12.  Effect  of  Le/B  ratio  on  pile  rotation  at  soil-cap  interface.  HP  14X73  pile  with  15  kip  shear  at  base  of  cap 


Figure  16.  Effect  of  Le/B  ratio  on  pile  rotation  at  soil-cap  interface,  HP  14X73  pile  in  stiff  clay,  k  =  50  tons/cu  ft 


HARD  CLAY,  k=  200 ton/ft3 


Ratio  of  Embedment  Length  to  Pile  Width,  Le/B  (m./in.) 


Figure  19.  Effect  of  Lc,  B  ratio  on  pile  rotation  at  soil-cap  interface,  various  pile  sections  in  medium  dense  satul 
with  10  kip  shear  at  base  of  cap. 
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constant  value  independent  of  Le/B  which  indicates 
the  fixed  end  condition  is  achieved. 

Figures  10  through  13  show  the  effect  of  soil  type 
and  soil  stiffness.  As  the  stiffness  of  soil  increases,  the 
rotation  of  pile  at  the  soil-cap  interface  decreases  for 
a  given  value  of  Lc/B.  Figures  I  2  and  13  show  that  the 
rotation  of  the  pile  is  less  for  clayey  soil  than  sandy 
soil  for  a  given  value  of  Lo/B. 


Figures  14  through  17  show  the  effect  of  shear 
force  at  the  interface  between  the  soil  and  the  cap.  As 
the  shear  force  increases  for  a  given  value  of  Le/B,  the 
rotation  of  the  pile  at  the  soil-cap  interface  increases.  • 


The  effect  of  pile  stiffness  is  depicted  in  Figures  18 
and  1‘).  In  both  figures  various  HP  pile  shapes  are 
examined.  These  two  figures  show  that  as  the  pile  silli¬ 
ness  increases  the  rotation  of  the  pile  at  the  soil-cap 
interface  decreases  for  a  given  value  of  I.e  'B. 


Table  8  lists  the  computed  response  for  a  typical  de¬ 
sign  situation.  In  this  table,  shear  at  the  base  of  the  cap 
was  used  as  an  input  (from  Table  2)  to  compute  the 
maximum  moment  and  cap  displacement.  The  analysis 
was  performed  for  12  in.,  24  in..  36  in.,  and  48  in.  em¬ 
bedment  length,  and  two  different  piles.  HP  14X73 
and  HP  14X117.  The  last  column  of  the  table  shows 
the  fixed  end  design  moments  obtained  from  Table  2 
in  order  to  compare  them  with  the  computed  maxi¬ 
mum  moments.  It  should  be  noted  that  the  computed 
moments  for  the  12  in.  embedment  length  ranged  from 
61  to  83  percent  of  the  computed  moments  for  the  48 
in.  embedment  length.  This  suggests  that  the  "rule  of 
thumb"  commonly  used  in  the  design  of  HP  piles,  i.e.. 
an  embedment  length  of  1  ft  or  less  into  the  concrete 
structure  is  considered  pinned,  may  be  unrealistic,  and 
that  past  pile  foundations  designed  using  I  ft  embed¬ 
ment  length  for  the  pinned  assumption  may  not  really 
be  pinned,  but  partially  fixed. 


In  the  parametric  analyses,  a  fixed-end  condition 
was  achieved  whenever  the  embedment  length  in  the 
pile  head  equaled  or  exceeded  twice  the  depth  of  the 
pile.  This  fixed-end  condition  occurs  when  the  rotation 
at  the  pile-cap  interface  reaches  a  maximum  value  inde¬ 
pendent  of  any  further  increase  in  the  embedment  of 
the  pile. 


axial  stress  distribution  A 

4  CAP-MEMBER  MODELS 


d) 


bending  stress  distribution 


Figure  20.  Stress  distribution  in  an  HP  pile  section. 


Introduction 

tTERL  analyzed  the  degree  of  fixity  of  the  pile  em¬ 
bedded  in  a  concrete  cap  from  a  purely  structural 
standpoint,  without  considering  the  underlying  soil  and 
the  interaction  of  the  pile  and  the  surrounding  soil. 
CERL  investigated  the  behavior  of  steel  HP  pile  sec¬ 
tions  with  several  lengths  of  embedment  in  various  con¬ 
crete  caps  using  two-dimensional  linear  and  nonlinear 
finite  element  models.  The  moments,  shears,  and  axial 
loads  that  would  normally  exist  in  the  HP  pile  section 
at  the  concrete  cap  face  were  assumed  to  produce  the 
stress  distributions  in  the  HP  pile  sections  that  are  pre¬ 
sented  in  Figure  20.  These  stress  distributions  were 
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used  to  determine  the  nodal  forces  to  apply  to  the 
appropriate  nodal  points  in  the  finite  element  models. 

Exploratory  Finite  Element  Models 

Before  constructing  the  finite  element  models  for 
the  typical  design  loading  conditions  and  the  para¬ 
metric  analyses,  four  exploratory  two-dimensional 
linear  models  of  a  concrete  cap  and  embedded  pile 
section  were  developed.  The  purpose  was  to  (a)  obtain 
a  general  understanding  of  the  fundamental  behavior  of 
the  concrete  cap  and  embedded  pile  section,  (b)  locate 
zones  of  stress  concentration,  and  (c)  investigate  the 
degree  of  detail  required  to  adequately  model  the  con¬ 
crete  cap  and  pile  section.  The  four  exploratory  models 
represented  an  HP  12X74  pile  embedded  a  depth  of 
2  ft  in  a  concrete  cap  with  a  nominal  width  of  36  in. 
(3D)  and  a  depth  of  48  in.  (4D).  where  the  symbol  D 
represents  the  depth  of  the  pile  section.  The  thickness 
normal  to  the  two-dimensional  plane  was  12.2  in.,  i.e.. 
the  width  of  the  flange  of  a  HP  12X74  pile.  Quadri¬ 
lateral  nonconforming  elements  were  used  to  represent 
the  embedded  pile  section  and  the  concrete  cap.  The 
aspect  ratio  for  the  elements  was  chosen  to  ensure  that 
the  resultant  aspect  ratios  were  not  greater  than  5:1 
and  that  the  majority  of  the  elements  had  ratios  no 
greater  than  2:1.  The  finite  element  meshes  for  the 
four  models  ranged  from  coarse  to  fine  in  an  effort  to 
determine  an  optimum  mesh.  The  optimum  mesh 
would  be  the  coarsest  mesh  that  would  yield  rea¬ 
sonable  results  as  compared  with  the  results  from  the 
finest  finite  mesh.  The  four  finite  element  models  are 
shown  in  Figures  21.  The  total  number  of  elements  in 
each  of  the  four  models  ranging  from  coarsest  to  finest 
was  12,  48.  350.  and  800.  respectively.  In  the  two  finer 
mesh  models,  the  thickness  of  the  flange  of  the  HP 
12X74  pile  was  modeled  as  a  single  row  of  elements. 
The  boundary  conditions  at  the  right  and  left  bound¬ 
aries  of  the  model  were  restrained  in  both  the  vertical 
and  horizontal  directions,  simulating  a  condition  of 
fixity  at  these  boundaries. 

The  behavior  of  both  the  concrete  cap  and  the  steel 
HP  pile  sections  were  assumed  to  be  governed  by  linear 
elasticity,  although  it  was  recognized  that  the  behavior 
of  concrete  is  nonlinear.  A  perfect  bond  was  assumed 
to  exist  between  the  pile  section  and  the  concrete  cap, 
i.e.,  displacements  at  the  interface  between  the  con¬ 
crete  cap  and  the  pile  were  identical  in  each  material. 
The  cross  sectional  area  of  the  pile  section  was  trans¬ 
formed  to  an  equivalent  concrete  thickness  using  the 
ratio  of  the  modulus  of  elasticity  for  steel  and  con¬ 
crete.  The  nodal  forces  applied  to  the  nodal  points  for 
the  embedded  pile  represented  a  constant  horizontal 


shear  force  of  I  kip  and  an  applied  moment  of  ’’20 
in. -kips. 

The  four  two-dimensional  finite  element  models 
were  analyzed  using  the  ANSYS  finite  element  com¬ 
puter  program.  Displacement  fields  for  these  models 
appeared  to  coincide  very  closely,  suggesting  that  the 
model  with  the  coarsest  mesh  would  be  appropriate. 
The  stress  distribution,  however,  did  not  substantiate 
the  same  conclusion.  For  a  nonconforming  quadri¬ 
lateral  element,  the  stresses  are  constant  within  the 
boundaries  of  the  element.  Therefore,  the  more  re¬ 
fined  the  mesh,  the  more  detailed  the  stress  distribu¬ 
tion.  When  the  mesh  becomes  very  fine  the  distribution 
of  stresses  reflects  ''peak"  stresses,  and  when  the  mesh 
becomes  coarse  the  stress  distribution  reflects  "aver¬ 
age"  stresses.  Based  on  a  comparison  of  the  stress  dis¬ 
tribution  for  the  four  models,  it  was  concluded  that 
the  second  model  (see  Figure  21b)  was  adequate  lor 
the  present  study. 

In  addition,  a  series  of  analyses  was  performed  w  ith 
the  mesh  shown  in  Figure  21c  to  gain  a  basic  under¬ 
standing  of  the  effect  of  embedment  depth  on  the 
magnitude  of  the  rotation  of  the  pile  section  at  the 
face  of  the  concrete  cap.  The  depths  of  embedment  of 
the  pile  were  0.5D,  D,  2D.  and  4D  or  6  in..  12  in..  24 
in.,  and  48  in.,  respectively.  The  nodal  forces  applied 
to  the  mode)  represented  a  horizontal  shear  force  of 
1.0  kip  and  applied  moments  of  240.  480.  and  720 
in.-kip.  The  rotation  of  the  pile  section  at  the  face  of 
the  concrete  cap  is  listed  in  Table  9  as  a  function  of  the 
applied  loads  and  the  depth  of  embedment.  The  lateral 
displacement  curves  for  the  centerline  of  the  HP 
12X74  pile  under  the  three  loading  conditions  with 
embedment  depths  of  0.5D.  D.  2D.  and  4D  are  pre¬ 
sented  in  Figure  22. 

As  the  depth  of  pile  embedment  is  increased  for  a 
given  pile  section  and  loading  conditions,  the  angular 
rotation  at  the  concrete  cap  face  should  gradually  de¬ 
crease  until  full  Fixity  is  achieved.  This  condition  will 
occur  when  a  limiting  value  of  angular  rotation  is 
reached  for  which  there  is  no  further  reduction  in 
angular  rotation  associated  with  additional  increases  in 
the  depth  of  pile  embedment.  For  the  rotational  results 
in  Table  9.  the  largest  decrease  in  angular  rotation  at 
the  concrete  cap  face  occurred  when  the  pile  embed¬ 
ment  depth  was  increased  from  0.5  to  ID  and  ranged 
from  12  to  25  percent  depending  on  the  loading  con¬ 
dition.  When  the  pile  embedment  depth  was  increased 
from  1  to  2D.  only  about  a  5  percent  decrease  in  the 
angular  rotation  was  observed  for  all  the  loading  con¬ 
ditions.  These  results  indicate  that  for  each  of  the  three 
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Figure  21.  Exploratory  finite  element  models. 


loading  conditions,  the  angular  rotation  has  reached 
some  limiting  value  at  an  embedment  depth  of  ap¬ 
proximately  2D  or  greater.  For  pile  embedment  depths 
less  than  approximately  2D,  only  partial  fixity  was 
achieved. 

Because  concrete  cracks  at  low  tensile  stress  levels 
and  degradation  occurs  under  cyclic  loading,  it  was 
decided  that  the  current  models  did  not  represent  the 
load  transfer  mechanism  at  the  interface  between  the 
concrete  cap  and  the  pile  section.  Therefore,  an  inter¬ 
face  element  was  incorporated  into  the  model  pre¬ 
sented  in  Figure  21b  before  proceeding  with  the 
additional  investigation  discussed  below. 

HP  14X117  Typical  Design  Loading  Analyses 

After  the  analyses  of  the  exploratory  finite  element 
models  were  completed,  the  St.  Louis  District  asked 
CERL  to  model  an  HP  14X117  pile  embedded  in  a 
concrete  cap  with  2-ft  and  4- ft  embedments  (nominal 
2D  and  4D  embedments).  These  models  were  to  be 
subjected  to  the  typical  loading  from  an  actual  design 
condition  as  presented  in  Table  2.  The  knowledge 
gained  from  the  exploratory  analyses,  in  particular  the 
need  to  incorporate  an  interface  clement,  was  factored 
into  the  development  of  the  finite  element  models. 
Specific  details  regarding  the  models  are  summarized 
below: 

1.  Pile-Quadrilateral,  nonconforming  elements 
were  used  to  model  the  web  of  the  pile,  plane  bar  ele¬ 
ments,  STIF1,  were  used  to  model  the  flanges  of  the 
pile.  The  thickness  of  the  finite  elements  was  equi¬ 
valent  to  the  thickness  of  web  of  the  pile  -0.805  in. 
To  simulate  the  flanges,  plane  bar  elements  with  an 
area  equal  to  the  area  of  the  flange-11.982  sq  in. 
were  connected  to  the  vertical  nodes  on  the  extreme 
right  and  left  sides  of  the  finite  elements  representing 
the  web  (Figure  23).  Linear  elastic  properties  were  as¬ 
sumed  for  the  steel.  The  modulus  of  elasticity  and 
Poisson's  ratio  were  assumed  to  be  29.6  X  106  lb/sq 
in.  and  0.3,  respectively. 

2.  Concrete  cap -Quadrilateral  nonconforming  ele¬ 
ments  were  also  used  to  model  the  concrete  cap.  The 
thickness  of  the  elements  outside  the  region  of  the  pile 
was  taken  as  the  width  of  the  flange  of  the  HP  14X1 17 
pile  a  conservative  assumption  for  this  analysis.  The 
width  of  the  flange  was  14.885  in.  The  thickness  of  the 
elements  representing  the  concrete  between  the  flanges 
of  the  HP  pile  section  was  reduced  by  the  thickness  of 
the  web  of  the  HP  pile  section  0.805  in.  In  essence, 
the  elements  representing  the  concrete  between  the 
flanges  and  the  elements  representing  the  web  of  the 


IIP  pile  section  overlaid  each  other  and  had  common 
nodal  points  at  the  corners  o|  the  elements,  the  con- 
ciele  between  the  llaugcs  ot  the  pile  would  delorm  in 
the  same  manner  as  the  IIP  pile  section.  I.meai  elastic 
properties  were  also  assumed  loi  the  conciete.  The 
modulus  ol  elasticity  and  Poisson  s  ratio  were  assumed 
to  be  3.0  X  I  <>*  lb  sq  in.  and  0  15.  respect ively  . 

3.  Pile  and  concrete  cap  interface  To  simulate  the 
load  transfer  mechanisms  between  the  external  sur¬ 
faces  of  the  flanges  of  the  HP  pile  section  and  the  con¬ 
crete  cap.  as  well  as  with  the  end  of  the  HP  pile  section 
and  the  concrete  cap,  a  two-dimensional  interface  ele¬ 
ment  (see  Figure  24)  was  introduced  into  the  model. 
This  element  is  designated  ST  IF  1  2  in  the  ANSYS  pro¬ 
gram.  It  represents  two  plane  surfaces  which  may  main¬ 
tain  or  break  physical  contact  and  may  slide  relative  to 
each  other.  The  element  is  capable  of  supporting  com¬ 
pression  in  the  direction  normal  to  the  surfaces  and 
shear  (Coulomb  friction)  in  the  tangential  direction. 
The  basic  features  and  force  deflection  relationships  of 
the  interface  elements  are  illustrated  in  Figure  25. 

The  force  deflection  relationships  for  the  interface 
elements  can  be  separated  into  the  normal  and  tangen¬ 
tial  directions  as  shown  in  Figure  25.  The  element 
operates  bilinearly  and  requires  an  iterative  solution. 
The  stiffness  matrix  is  reformulated  after  each  itera¬ 
tion  for  static  convergence. 

The  element  status  at  the  beginning  of  the  first 
iteration  was  assumed  to  be  closed;  the  stiffness,  k.  was 
used  for  computing  both  gap  and  sliding  effects.  The 
element  status  at  the  beginning  of  each  succeeding 
iteration  was  based  on  the  status  at  the  end  of  the  pre¬ 
vious  iteration.  In  the  normal  direction,  if  the  normal 
force  (Fn)  was  negative,  the  interface  remained  in  con¬ 
tact  and  responded  as  a  linear  spring.  If  the  normal 
force  became  positive,  contact  was  broken  and  no 
force  was  transmitted  across  the  interface.  In  the  tan¬ 
gential  direction,  if  Fn  <0  and  the  absolute  value  of 
the  tangential  forces,  Fs,  was  less  than  or  equal  to  the 
interface  coefficient  of  friction  (p).  times  the  absolute 
value  of  fn  (p|FnQ,  the  interface  did  not  slide  and 
responded  as  a  linear  spring  in  the  tangential  direction. 
However,  if  Fn  <  0  and  Fs  >  pFn.  sliding  occurred. 
If  Fn  >  0  contact  was  broken  and  the  tangential  force 
was  reduced  to  zero.  Convergence  occurred  whenever 
the  status  of  the  bilinear  elements  remained  unchanged 
from  the  previous  iterations. 

The  only  material  property  required  lot  the  inter¬ 
lace  element  is  the  coefficient  of  friction,  p.  The  P*e 
stressed  Concrete  Institute  (PCI)  specifies  a  coelticient 
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Figure  25.  Force-deflection  relationship  for  the  interface  element. 


of  friction  of  0.4  for  a  concrete  and  steel  interface.'0 
The  American  Concrete  Institute  (AC1)  specifies  a 
coefficient  of  friction  of  0.7  when  surfaces  of  rolled 
steel  and  concrete  are  interacting.11  The  AC  I  also 
specifies  a  coefficient  of  friction  of  1.4  for  concrete 
placed  monolithically,  which  is  an  upper  limit  for  the 
coefficient  of  friction.  A  coefficient  of  friction  of  0.4 
was  selected  for  the  initial  analyses,  because  it  repre¬ 
sented  a  lower  limit. 

4.  Boundary  conditions  The  same  boundary  con¬ 
ditions  were  used  for  the  HP  14X1 1  7  pile  model  as  for 
the  HP  1  2X74  pile  model. 

5.  Loading  The  loading  conditions  for  the  HP 
14X117  pile  section  are  shown  in  Table  2.  The  intro¬ 
duction  of  the  interface  element  required  that  the 
loads  be  gradually  applied  to  the  model  to  ensure  con¬ 
vergence  of  the  solution.  Consequently,  the  shear,  a\i. ! 
load,  and  moment  for  each  of  the  loading  conditions 
were  gradually  applied  to  the  model  in  a  series  of  load 
steps  so  that  large  corrections  in  succeeding  iterations 
would  not  be  required.  After  each  load  step  the  solu¬ 
tion  was  allowed  to  iterate  until  convergence  was  ob¬ 
tained  as  represented  by  the  unchanged  status  of  the 
bilinear  elements  from  the  previous  iteration. 

In  the  analysis  of  the  HP  14X1 17  pile  section  with 
2D  and  4D  embedments  and  the  subsequent  parametric 
analysis  (discussed  below),  the  effects  of  various 
parameters  were  analyzed.  Table  10  is  an  index  correla¬ 
ting  the  analysis  case  number  and  the  specific  model 
parameters  for  the  analyses  summarized  in  this  report. 
Table  10  should  be  used  as  a  reference  when  reviewing 
and  comparing  some  of  the  figures  which  follow. 

A  typical  displacement  field  and  stress  contour  plots 
for  the  HP  14X117  pile  under  normal  operations  with 
a  0. 1  g  earthquake  condition  and  an  embedment  depth 
of  2D  are  presented  in  Figures  26  through  32.  The 
displacement  field  (Figure  26)  indicates  that  although 
the  maximum  displacement  was  only  about  0.004  in., 
the  pile  section  and  concrete  cap  separated  at  the  right 
side  of  the  pile  section.  Consequently,  relative  displace¬ 
ments  occurred  between  the  concrete  cap  and  pile 
section  along  both  vertical  interfaces,  as  shown  by  the 
offset  boundaries  of  the  finite  elements  along  these 
interfaces. 


10 PCI  Manual  on  Design  of  Connections  lor  Precast.  Pre 
stressed  Concrete  (Prcsiressed  Concrete  Institute.  1973) 
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Stress  contour  plots  for  ox.  ox .  txx.  oinax.  oIIljn . 
and  tmax  are  presented  in  Figures  27  through  32.  In 
each  normal  stress  contour  plot  (Figs.  27.  23.  30.  31  ). 
maximum  tensile  and  compressive  stress  levels  and 
their  points  of  occurrence  are  indicated  with  an  "X” 
and  “0,"  respectively.  In  the  shear  stress  contour  plots 
(Figs.  20,  32),  the  maximum  and  minimum  stress  levels 
and  their  points  of  occurrence  are  similarly  marked. 
The  stress  contour  interval  is  indicated  near  the  lower 
left  corner  of  the  plot.  A  dash  line  is  used  to  represent 
the  zero  stress  contour  I, ne.  Figures  27  and  2k  indicate 
that  the  maximum  ox  and  ux  stresses  were  30S  Ib/sq 
in.  and  371  lb/sq  in.,  respectively.  These  stresses 
occurred  at  the  concrete  cap  face  at  the  interface  be¬ 
tween  the  concrete  cap  and  pile  on  its  left  side.  Along 
the  interface  between  the  concrete  cap  and  the  left  side 
of  the  pile,  the  ox  stresses  gradually  decreased  to  a 
stress  level  approaching  zero  near  the  end  of  the  em¬ 
bedded  pile.  Along  this  same  interface,  the  stresses 
(see  Figure  2K)  decreased  much  more  slowly  because 
additional  load  was  transferred  from  the  pile  section  to 
the  concrete  cap  through  the  interface  elements. 

Figures  30  through  32  indicate  that  the  maximum 
and  minimum  principal  stresses  occurred  at  the  left 
side  of  the  pile  section  at  the  concrete  cap  face  and 
that  the  maximum  shear  force  occurred  along  the  left 
boundary  of  the  concrete  cap.  The  maximum  tensile 
stress  was  106  Ib/sq  in.  and  the  maximum  shear  stress 
was  234  Ib/sq  in.  The  latter  is  slightly  greater  than 
i.S\/Q  which  is  the  maximum  allowed  by  the  AC' I 
code  lor  members  subjected  to  shear  and  flexure 
only.12  Since  it  is  standard  construction  practice  to 
provide  concrete  steel  reinforcement  only  in  the  region 
above  the  top  of  the  embeuded  pile,  consideration 
should  be  given  to  placement  of  a  nominal  percentage 
of  concrete  steel  reinforcement  in  the  region  of  pile 
embedment  to  increase  the  structural  integrity  of  the 
concrete. 

The  lateral  displacements  of  the  centerline  of  the 
pile  sections  with  4D  and  2D  embedment  in  the  con¬ 
crete  cap  (Case  1  and  <^ase  2.  respectively)  are  shown 
in  Figure  33.  The  displacement  curve  for  2D  embed¬ 
ment  nearly  coincided  with  the  displacement  curve  lor 
4D  embedment  between  0  and  2D.  In  addition,  the 
displacement  curve  for  4D  embedment  experienced  a 
change  in  slope  at  a  distance  slightly  gieater  than  2D 
from  the  concrete  cap  face,  which  indicates  that  the 
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Table  10 

Index  Correlating  Analysis  Case  and  Model  Parameters 
Model  Parameters 


Fqgure  26.  Typical  displacement  field  for  HP  14X117  pile 
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Figure  31.  Typical  Stress  Contour  Plot  for 


Figure  32.  Typical  Stress  Contour  Plot  for 


pile  section  has  achieved  a  condition  of  full  fixity.  The 
maximum  lateral  displacement  at  the  face  of  concrete 
cap  for  the  pile  section  with  4D  embedment  was 
slightly  greater  than  the  pile  section  with  2D  em¬ 
bedment.  Considering  the  nonlinear  aspects  of  the 
concrete  pile  system,  this  difference  is  not  considered 
significant. 

The  rotations  of  the  pile  section  at  the  face  of  the 
concrete  cap,  </>max,  for  these  same  conditions  of  2D 
and  4D  embedment  were  computed  to  be  3.22  X  10'4 
rad  and  3.20  X  I04  rad,  respectively.  A  difference  of 
about  2  percent  between  the  maximum  rotations  for 
the  2D  and  4D  embedments  was  not  considered  signifi¬ 
cant.  Also,  the  magnitudes  of  these  rotations  were 
quite  similar  to  those  obtained  in  the  analysis  of  the 
exploratory  finite  element  models  for  the  case  of  2D 
embedment  and  a  bending  moment  of  720  in. -kip. 

The  compression  and  shear  forces  developed  in  the 
interface  elements  in  the  models  are  shown  in  Figures 
34  and  35  for  the  2D  and  4D  embedments,  respec¬ 
tively,  together  with  the  nodal  forces  simulating  the 
applied  loads.  The  forces  transmitted  by  the  first  two 
interface  elements  on  the  lower  left  side  of  the  pile  sec¬ 
tion  were  nearly  identical  for  the  2D  and  4D  embed¬ 
ments.  Forces  in  the  other  interface  elements  along  the 
left  side  of  the  pile  section  tended  to  degenerate  to 
zero,  thus  indicating  that  contact  between  the  pile  sec¬ 
tion  and  the  concrete  cap  had  been  broken.  This  sug¬ 
gests  that  the  pile  section  attempts  to  pivot  about  the 
edge  of  the  concrete  cap  and  thereby  reduce  the  forces 
acting  on  the  upper  portion  of  the  left  side  of  the  pile 
section;  this  increases  the  forces  acting  on  the  upper 
portion  of  the  right  side  of  the  pile. 

Due  to  the  nonlinear  aspects  of  the  interface  ele¬ 
ments  and  the  solution  convergence  criteria,  the  forces 
acting  at  the  interface  elements  do  not  fully  satisfy 
statics.  However,  the  estimated  error  was  less  than  4 
percent. 

HP  14X117  Parametric  Analysis 

A  parametric  analysis  of  the  HP  14X117  pile  em¬ 
bedded  2D  in  a  concrete  cap  was  performed  to  deter¬ 
mine  if  other  parameters  would  significantly  affect  the 
degree  of  fixity  of  the  pile  section.  The  loadings  ap¬ 
plied  to  the  pile  were  those  associated  with  normal 
operations  with  0.1  g  earthquake  conditions.  The  basic 
parameters  for  the  models  were  the  same  as  those  used 
in  previous  analysis  and  are  summarized  in  Table  11. 
The  following  parameters  were  investigated. 


1 .  Coefficient  of  friction  between  the  concrete  cap 
and  the  pile  section 

2.  Effective  width  of  the  concrete  cap 

3.  Pile  spacing 

4.  Concrete  cap  thickness. 

Different  finite  element  models  were  developed  and 
analyzed  based  in  part  on  the  range  of  parameters  pre¬ 
sented  in  Table  1 .  Table  10  is  an  index  which  correlates 
the  model  parameters  with  a  specific  analysis  case. 

The  following  is  a  summary  of  the  results  from  the 
parametric  analysis: 

1 .  Coefficient  of  friction  between  the  concrete  cap 
and  the  pile  section.  A  degree  of  adhesion  exists  be¬ 
tween  the  concrete  cap  and  the  pile  section.  This  ad¬ 
hesion  occurs  at  the  interface  between  the  concrete  cap 
and  the  exterior  surface  of  the  flanges  of  the  pile  sec¬ 
tion,  and  at  the  concrete  cap  and  the  end  of  the  pile. 
Consequently,  three  different  values  for  the  coefficient 
of  friction  were  investigated:  0.4,  0.7,  and  1.4.  Three 
analyses  were  performed  and  identified  as  Case  2  (p  = 
0.4),  Case  3  (/t  =  0.7),  and  Case  4  (p  =  1 .4). 

The  lateral  displacement  curves  for  the  centerline  of 
the  embedded  pile  section  for  the  three  analysis  are 
shown  in  Figure  36.  The  shapes  of  the  three  displace¬ 
ment  curves  arc  basically  similar.  However,  as  the 
coefficient  of  friction  was  increased  from  0.4  to  1 .4. 
the  maximum  displacement  of  the  centerline  of  the 
pile  section  at  the  concrete  cap  face  decreased  from 
2.4l>  X  10  3  in.  to  1 .62  X  10  3  in.,  or  about  35  per¬ 
cent.  The  maximum  rotation  decreased  from  3.22  X 
ICC4  rad  to  2.23  X  10  4  rad.  or  about  30  percent. 
These  effects  were  due  to  the  transfer  of  forces  in  the 
pile  section  to  the  concrete  cap  due  to  the  increase  in 
the  coefficient  of  friction.  Comparison  of  the  interface 
element  forces  in  Figures  34  (Case  2).  37  (Case  3).  and 
38  (Case  4)  clearly  indicated  this  to  be  the  load  trans¬ 
fer  mechanism. 

Stresses  in  the  concrete  cap  were  also  compared  at 
the  sections  indicated  in  Figure  39.  As  a  result  of  the 
nonconforming  quadrilateral  element  used  to  repre¬ 
sent  the  concrete  cap.  the  stress  obtained  from  the 
analyses  was  average  stress  at  the  center  of  the  element 
The  distribution  of  forces  in  the  horizontal  direction 
(fx)  and  the  vertical  direction  (fy )  were  calculated 
using  the  normal  stress,  ox ;  the  shear  stress.  tvv  .  and 
the  thickness  of  the  concrete  cap,  14.^  in.  These  forces 
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Lateral  Displacement,  I0~4in. 

Figure  33.  Lateral  displacement  curves  for  centerline  of  HP  14X1 17  pile  with  2D  and  4D  embedments,  p  =  0.4. 
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Figure  34.  Interface  element  forces  for  HP  14X1 1  7  pile  with  2D  embedment,  p  =  0.4. 
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Figure  35.  Interface  element  forces  for  HP  14X117 
pile  with  4D  embedment,  /a  =  0.4. 


were  plotted  in  Figures  40  through  43  In  general,  the 
force  distributions  are  similar  to  each  other,  which 
indicates  that  the  concrete  cap  stresses  were  not  signifi¬ 
cantly  influenced  by  the  coefficient  of  friction,  at  least 
for  the  range  of  values  that  were  investigated. 

2.  Effective  width  of  the  concrete  cap.  The  effec¬ 
tive  width  of  the  concrete  cap  was  increased  from 
143)  in.  to  44.7  in.  to  determine  if  the  effective  width 
of  the  concrete  cap  significantly  affected  the  behavior 
of  the  cap-member  system.  To  model  this  condition,  the 


thickness  of  the  elements  representing  the  conciete  s.ip 
outside  the  region  of  the  pile  was  increased  to  44  "  in 
while  the  thickness  of  the  elements  represennne  the 
conciete  between  the  flanges  of  the  pile  section  was  m 
creased  to  43.9  in.  lhe  lallei  number  repiesenis  ihc 
effective  width  of  44. 7  in.lesslhelhickness.il  lhe  web 
of  the  pile  section.  Analysis  Case  5  was  peiloimed  and 
compared  with  the  results  from  Case  2.  The  Ijlcial 
displacement  curves  for  the  centerline  of  the  embed¬ 
ded  pile  section  with  an  effective  width  ot  concrete 
cap  of  14.9  in  and  44.7  in.  are  presented  in  Figuie  44 
The  maximum  displacement  at  the  centerline  ot  tin 
pile  section  at  the  concrete  cap  face  was  reduced  from 
2.49  X  10  4  in.  to  1.16  X  10  4  in.,  or  about  50  per 
cent.  The  maximum  rotation  at  the  concrete  cap  face 
was  reduced  from  3.22  X  10  4  rad.  to  I  "6  /  10 
rad  or  about  45  percent,  with  the  increase  in  the  effec¬ 
tive  width  of  the  concrete  cap.  The  shape  of  the  Case 
5  displacement  curve  was  different  in  that  the  slope  > -t 
the  curve  changed  sign  at  an  embedment  depth  ot 
about  I.75D.  Likewise  the  distributed  forces  shown  in 
Figures  40  through  43  for  Case  2  and  Case  5  do  not 
indicate  significant  differences. 

3.  Pile  spacings.  Pile  spacings  varying  between  3D 
and  5D  were  specified  in  Table  1.  Analy  sis  Case  6  was 
performed  to  reflect  the  specified  variation  in  the  pile 
spacing  parameter  and  compared  with  the  results  trom 
Case  2. 

The  lateral  displacement  curves  for  the  embedded 
pile  section  with  pile  spacings  of  3D  (Case  2)  and  5D 
(Case  6)  are  presented  in  Figure  45.  The  displacement 
curves  are  essentially  identical.  When  the  pile  spacing 
was  increased,  the  maximum  lateral  displacement  in¬ 
creased  from  2.4‘)  X  10  4  in.  to  2.76  X  10  4  in.,  oi 
about  10.5  percent,  and  the  maximum  rotation  in¬ 
creased  from  3.22  X  10  4  rad  to  3.31  X  10  4  rad.  or 
about  2.6  percent.  Comparison  of  the  distributed  force 
plots  for  fx  and  t\  in  Figures  40  thtough  43  ami 
Figures  46  through  4l>  revealed  that  there  were  only 
minor  differences  between  the  two  analyses 

4.  Concrete  Cap  Thickness.  A  series  of  analyses  was 
performed  in  which  the  pile  section  was  embedded  2D 
in  the  concrete  cap  and  the  thickness  ot  concrete  covei 
above  the  pile  section  was  varied  from  41)  to  ~D 
(represented  as  Case  7.  Case  8.  and  Case  9).  Compari¬ 
sons  of  the  distributed  forces  calculated  liom  the 
results  of  these  analyses  are  presented  in  Figures  4o 
through  4».  respectively.  These  plots  show  that  the 
thickness  of  concrete  cover  above  the  pile  section  does 
not  alter  the  distribution  of  forces  in  the  legion  ot  the 
depth  of  embedment  for  the  pile  section 
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Table  1 1 

Basic  Parameters  for  HP  14X1 17 
Pile  Parametric  Analysis 


Parameter 

Value 

Member  Spacing 

31) 

Member  [imbedment 

2D 

Cap  Thickness 

4D 

Cap  Effective  Width 

14.9  in. 

Member  Configuration 

HP  14X117 

Modulus  of  l Elasticity 

Concrete 

3.0  x  10‘  Ib/sq  in. 

Steel 

29.6  x  10‘  Ib/sq  in 

Poisson  Ratio 

Concrete 

0.15 

Steel 

0.3 

Coefficient  of  friction 

0.4 
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Case  4 


Figure  37.  Interlace  element  forces  for  HP  14X1 17  pile  with  2D  embedment.  »  -  0.7. 
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Figure  39.  Sections  in  concrete  cap  used  for  normal 
and  shear  force  distribution  comparisons. 


5  CONCLUSIONS 


Based  on  the  results  from  the  analytical  studies  by 
UIUC  and  CERL  on  the  cap-member-soil  system  and 
the  cap-member  system,  a  pile  embedment  length  equal 
to  or  greater  than  twice  the  pile  depth  or  diameter  is 
required  to  develop  a  condition  approximating  full 
fixity  for  a  pile  embedded  in  the  base  of  the  structure. 
Moreover,  the  "rule  of  thumb”  commonly  used  in  the 
design  of  HP  piles  that  I  ft  embedment  length  into  the 


concrete  structure  is  considered  pinned  may  be  un¬ 
realistic.  The  computed  response  of  die  cap-nicinber- 
soil  model  for  the  typical  design  •>:  aiions  indicated  a 
1  ft  embedment  length  would  develop  (■>!  to  83  percent 
of  the  moments  for  a  4  ft  embedment  length.  There¬ 
fore  existing  pile  foundations  designed  using  I  ft  em¬ 
bedment  length  into  the  concrete  structure  for  the 
pinned  assumption  may  not  really  be  pinned,  but  par¬ 
tially  fixed. 
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Figure  40.  Norma!  force  distribution  at  Section  A-A  for  cap  thickness  of  4D, 
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